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I Abstract We study the asymptotic behavior of distribution densities arising in stock price models with 

^ ' stochastic volatility. The main objects of our interest in the present paper are the density of time averages 

^ I of the squared volatility process and the density of the stock price process in the Stein-Stein and the Heston 

model. We find explicit formulas for leading terms in asymptotic expansions of these densities and give 
error estimates. As an application of our results, sharp asymptotic formulas for the implied volatility in the 
Stein-Stein and the Heston model are obtained. 
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O^; 1 Introduction 

In [11] and [12], we found sharp asymptotic formulas for the distribution density of the stock price process 
^ , in the Hull- White model. These formulas were used in [13] to characterize the asymptotic behavior of the 

^SJ ' implied volatility in the Hull- White model. The present paper is a continuation of [11], [12], and [13]. It 

On , concerns the asymptotic behavior of various distribution densities arising in the Stein-Stein and the Heston 

' model. 

The stochastic differential equations characterizing the stock price process Xt and the volatility process 
' Yt in the Stein-Stein model have the following form: 

o- 

^ ■ / rfAt - ^iXtdt +\Yt\ XtdWt , . 

dYt^q{Tn-Yt)dt + adZt, ^ ' 



O 



where Wt and Zt are standard Brownian motions. The initial condition for Xt is denoted by xq and for Yt 
by yo. The stochastic volatility model in ^ was introduced and studied in [T^]. In this model, the absolute 
■ value of an Ornstein-Uhlenbeck process plays the role of the volatility of the stock. We will assume that that 

fi eR, q > 0, m> 0, and a > 0. 

The Heston model was introduced in [TJ]. In this model, the stock price process Xt and the volatility 
process Yt satisfy the following system of stochastic differential equations: 

dXt - fiXtdt + VYtXtdWt , . 

dYt ^ {a + bYt) dt + c^YtdZt , ^ ' 

where Wt and Zt are standard Brownian motions. We will assume that /i G M, a > 0, 6 < 0, and c > 0. The 
initial conditions for Xt and Yt are denoted by Xq and yo, respectively. The volatility equation in ^ can be 
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rewritten in the mean-reverting form. This gives 

r dXt ^ ^iXtdt + VYtXtdWt 

\ dYt^r{m- Yt) dt + c^YtdZt, 

where r = —b and m — The volatihty equation in Q is uniquely solvable in the strong sense, and 

the solution Yt is a positive stochastic process. This process is called a Cox-Ingersoll-Ross process (a CIR- 
process) . This process was studied in [6j . Interesting results concerning the Heston model were obtained in 

m- 

It will be assumed throughout the present paper that the models described by ^ and ^ are uncorrelated. 
This means that the Brownian motions Wt and Zt driving the stock price equation and the volatility equation 
in ^ and ^ are independent. In the analysis of the probability distribution of the stock price Xt, the 
mean-square averages of the volatility process over finite time intervals play an important role. For the 
Stein-Stein model, we set 

Here Yt satisfies the second equation in ([I]). For the Heston model, we have 

where Yt satisfies the second equation in It will be shown below that for every t > 0, the probability 
distribution of the random variable at defined by ([3]) for the model in ([T]) and by (|4]) for the model in ([2|) 
admits a distribution density mt (see Lemma 16.31 and Remark 18.11 below). The function mt is called the 
mixing distribution density. 

The distribution density of the stock price Xt will be denoted by D*. The existence of the density Dt 
follows from formula In this paper, we obtain sharp asymptotic formulas for the distribution density Dt 
in the case of the uncorrelated Stein-Stein and Heston models. Note that in [7] and [19], the behavior of Dt 
was studied for the Heston model and the Stein-Stein model, respectively, using a rough logarithmic scale in 
the asymptotic formulas. Moreover, no error estimates were given in these papers. The results established 
in the present paper are considerably sharper. We find explicit formulas for leading terms in asymptotic 
expansions of Dt in the Heston and the Stein-Stein model with error estimates. It would be interesting to 
obtain similar results for correlated models, since such models have more applications in finance. We hope 
that the methods employed in the present paper may be useful in the study of the correlated case. 

We will next quickly overview the structure of the present paper. In Section [21 the main results of the 
paper (theorems 12. 11 [2T2l 12.31 and l2.4p are formulated. They concern the asymptotic behavior of the mixing 
distribution density mt and the stock price distribution density Dt in the Heston model and the Stein-Stein 
model. Section [3| is devoted to applications of our main results. In this section, we obtain sharp asymptotic 
formulas for the implied volatility in the Heston and the Stein-Stein model. In section [4l we gather several 
known facts about CIR-processes and Bessel processes. We also formulate a theorem of Pitman and Yor 
concerning exponential functionals of Bessel processes. This theorem plays an important role in the present 
paper. In Section [5l we prove an asymptotic inversion theorem for the Laplace transform in a certain class 
of functions (see Theorem l5.2p . This theorem is useful in the study of the asymptotic behavior of the mixing 
distribution density, since it is often easier to find explicit formulas for the Laplace transform of this density 
than to characterize the density itself. In sections [6| - [U we prove theorems 12.11 - 12.41 and describe the 
constants appearing in these theorems. 

2 Asymptotic formulas for distribution densities 

The next four theorems are the main results of the present paper. The first two of them provide explicit 
formulas for the leading term in the asymptotic expansion of the distribution density Dt of the stock price 
Xt in the Stein-Stein model and the Heston model, while the other two concern the asymptotic behavior of 
the mixing distribution density nit in these models. 
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Theorem 2.1 Let Dt he the stock price distribution density in model (QP with q > 0, m > 0, and a > 0. 
Then there exist positive constants Bi, B2, and B3 such that 



Dt (xot 



Bix 



-B3 B2^/TogX 



(logx) 



1 + (log X 



(5) 



Theorem 2.2 Let Dt be the stock price distribution density in model (0j with a > 0, b < 0, and c > 0. 
Then there exist positive constants Ai, A2, and such that 



Dt {xoe^''x) = A^. 



/loR 



(6) 



Theorem 2.3 Suppose that q > 0, m > 0, and a > in model flj). Then there exist positive constants E 
F , and G such that 

mt{y) = Ee-^y\^y (l + O (y- 



y 



(7) 



Theorem 2.4 Suppose that a > 0, b < 0, and c > in model Then there exist positive constants A, B, 
and C such that 

00. (8) 



mt{y) = Ae-'^^'e^^y-^+f^ (l + O {y'^')) V 



It is known that for uncorrelated stochastic volatility models, the asymptotic behavior of the stock 
price distribution density Dt near zero is determined by the behavior of Dt near infinity. Indeed, we have 
Dt (xoe'^'x"^) — x^Dt (xoe'^'a;), a; > 0. This equality can be derived from the formula 



Dt 



(a;oe^*x) - 



1 



y ^mt{y)exp 



1,2 ty'^ 



dy 



(9) 



(see, e.g.. Section 4 in ^TT). It follows from Theorem 12.21 that the stock price distribution density Dt{x) in 
the Heston model behaves at infinity roughly as the function x~^^ and at zero as the function x^^~^ . For 
the Stein-Stein model. Theorem 12 . II implies that Dt{x) behaves at infinity roughly as the function x^^-^ and 



at zero as the function x^^~^ . In [19], the latter fact was established (in part, only heuristically) in the 
logarithmic scale (the same scale was used in [7]). More precisely, the following definition of the asymptotic 
equivalence was used in [7] and [19]: Two functions F{x) and G{x) are called asymptotically equivalent as 



X ^ 00 (or as X 0) if ^1 as x ^ 00 (or x — > 0). It is not hard to see that the constant 7 

in [19|, wich characterizes the decay of the stock price density Dt near infinity (see formula (20) in 19]), 
coincides with the constant i?3 in ([5]) (see lemmas 16.61 and 17.31 for the description of the constant B^). 
However, there is an error in formula (21) in [19 . It is stated that Dt{x) behaves near zero as the function 
2.-1+7^ By Theorem 12.11 the correct power function that characterizes the behavior of the density Dt near 



zero is the function 



-3+7 



The values of the constants A3 and B^ appearing in theorems 12.11 and 
formulas: 

A,^A,itAc) = l + ^^^^, C = CitAc)^^ 
3 VSGTt ^ , t 



and 



B3 = B3{t,b,c) 



are given by the following 
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2,-^2 

q +-^r^t 



(see lemmas 16.61 17.31 18.21 and 19.11 below) . Here denotes the smallest positive root of the entire function 
z I—!- zcosz + ssinz. The zeroes of this function are studied in Section [S| Note that A^ does not depend on 
a, while B3 does not depend on m. It is clear that A3 > 2 and B3 > 2. In [TT] and [H], we studied the 
tail behavior of the stock price distribution density Dt in the Hull- White model. It was established that the 
function Dt(x) behaves at infinitiy like the function x~^ on the power function scale. This is an extremely 
slow behavior. No uncorrelated stochastic volatility model has the function Dt{x) decaying like x"^"*"*^, e > 
(see [E]). The tail of the stock price distribution in the Hull- White model is "fatter" than the corresponding 
tail in the Heston and the Stein-Stein model. 
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3 Applications. Asymptotic behavior of the impUed volatiUty 



The implied volatility in an option pricing model is the volatility in the Black-Scholes model such that the 
corresponding Black-Scholes price of the option is equal to its price in the model under consideration. In 
the present paper, we will only consider European call options, and the implied volatility will be studied as 
a function of the strike price K. Let us denote by Vq the pricing function for the European call option in the 
Heston (or the Stein-Stein) model. Then the impled volatility I{K) satisfies Cbs{K,I{K)) = Vo{K), where 
Cbs stands for the Black-Scholes pricing function. 

In il3j, we studied the implied volatility in the Hull- White model. In the present paper, we characterize 
the asymptotic behavior of the implied volatility in the Stein-Stein and the Heston model. Our asymptotic 
formulas are sharper than the formulas which can be obtained form more general results due to Lee, Benaim, 
and Friz (see [15], [2], [2]). We will first consider the model in ([2]). It will be assumed below that the market 
price of volatility risk 7 is equal to zero (see, e.g., [9] for the definition of 7). Then, under the corresponding 
martingale measure P* , the system of stochastic differential equations in ([2]) can be rewritten in the following 
form: 

dXt = rXtdt + YtXtdW* 
dYt = {a + bYt) dt + cy/TtdZ; , 

where and are independent standard one-dimensional Brownian motions, and r > is a constant 
interest rate. This means that under the measure P* the new system also describes a Heston model. This is 
the reason why we assumed that 7 = 0. 

Let us consider a European call option associated with the stock price model in (fTO)) . The price of such 
an option at f = is given by the following formula: 

Vo{K)^E*[e-'^^{XT-K)^], (11) 

where T is the expiration date and K is the strike price. Put Dt{x) = Dt [x] r, v, ^, xq, yo). Then we have 



Vo{K) = e'''^ xDT{x)dx - e-"' K DT{x)dx. (12) 
Jk Jk 

The implied volatility is often considered as a function of the log-strike fc, which in our case is related to 
the strike price K by the formula k — log ^^^r'l - ■ In terms of fc, the implied volatility is defined as follows: 

I{k) ~ I{K), —00 < k < 00, < K < 00. For uncorrelated stochastic volatlity models, the behavior of 
the implied volatility near zero is completely determined by how it behaves near infinity. More precisely, it 
is known that the implied volatility is symmetric in the following sense: / (^{xoe'-'^y K-^^ =1{K) for all 

K > Q (see, e.g., [12]). It is clear that the previous equality can be formulated in terms of the log-strike k as 
follows: I{k) — I{~k), —00 < k < 00. 

Important results concerning the behavior of the implied volatility I{k) in a general case were obtained 
by Lee (see [15)). He characterized the asymptotic behavior of the implied volatlity for large strikes in terms 
of the moments of the stock price process. In [2], the asymptotic behavior of the implied volatility was linked 
to the tail behavior of the distribution of the stock price process (see also [3]). The reader interested in more 
aspects of the asymptotic behavior of the implied volatility can consult Chapter 5 in [9]. 

The following theorem will be established below: 



Theorem 3.1 For the Heston model 



as fc —> 00, where 



i{k)^p,kUp,+Ps'^+o(tm 

fc2 \ k^ J 

/3i = ^(v/i^-v/ir^), 



(13) 



2 



'2 



4 



^/2^(4 cO(vA^ 

xeasir; 

For the Stein-Stein model 



and ip can be any positive increasing function on (0, oo) such that lim ipik) = oo. 

k — ^oo 



as k —^ OO, where 



/(fc)=7ifcU72 + 0(^^) (14) 



71 = ^ (VS, - 1 - V^s - 2) , 



72 



/T 

B2 / 1 1 



and the function ijj is as above. 



4 CIR processes and Bessel processes 

The volatihty in model ^ is described by the square root of an CIR-process. It is known that CIR-processes 
are related to squared Bessel processes. In the present section, we gather several results concerning Bessel 
processes and CIR-processes. Let 6 > and x > 0, and consider the following stochastic differential equation: 
dTt = Sdt + 2\/j\dZt, Tq = X a.s. This equation has a unique nonnegative strong solution Tt, which is called 
the squared (S-dimensional Bessel process started at x. The following notation is often used for the squared 
Bessel process: Tt = BESQiit). We refer the reader to [H [3 [HI UHl [III HH] for more information on Bessel 
processes. 

The next lemma links Bessel processes and CIR-processes (see, e.g., TUj). 



Lemma 4.1 Let Yt be a CIR process satisfying the equation dYt = [a + bYt) dt + c^/YtdZt with Yq — x 

. 46 



-a.s., and put Tt = BESQf [t). Then Yt = e^*T (fj (l - e-^^] 



Remark 4.2 If 5 = 0, then we have Yt = BESQf (^i) . 

Pitman and Yor proved the following assertion (see pTl I18j ): 

Theorem 4.3 Let A > 0. Then 

\2 rt 



E 



exp|-y^ BESQiiu)du^ 



xX 

[cosh(At)] ^ exp ^ — — tanh(At) 



The next statement can be derived from Theorem 14.31 This statement can be found, e.g., in [4]. 



Theorem 4.4 Let a > 0, b < 0, c > 0, and let Yt be a CIR-process in Lemma \4.1\ such that Yq = yo a.s. 
Then for every rj > ^, 



exp 



abt 

„2 



f 277 \^ ^ \ byo (477^ - 1) sinh(bt77) \ 

\2r] coshlbtr]) — smhibtrj) J ^ | 2c'^ricosh{btr]) — c^ sm\i{btrj) j 



Theorem 14.41 is equivalent to the following assertion: 
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Theorem 4.5 Let a > 0, 6 < 0, c > 0, and let Yt be a CIR process in Lemma \4.1\ such that Yq — yo a.s. 
Then for every A > 0, 



E, 



Vo 



exp 



A / Y,ds 



= cxp 



exp • 




Vb^ + 2c2A 



2c2Acosh(kVfe2 + 2c2A) - 6sinh(iiV&2 + 2c2A) 



2?/oA sinh(itVfo2 + 2c2A) 



(15) 



\/PT2?Acosh(it%/62T2?A) - bsmh{^tVb^ + 2c^\) 

Theorem 14.51 characterizes the Laplace transform of the probability distribution of the random variable 
Ygds. In the next section, we will study the asymptotic behavior of the inverse Laplace transform in a 
class of functions which look like the function on the right-hand side of ([1 



5 An asymptotic inverse of the Laplace transform 

Let us assume that M is a function on (0, oo) whose Laplace transform is given by the following formula: 

e'^yM{y)dy = /(A), A > 0, 



where 

/(A) = A^iGi(A)''^G2(A)e-^(^) 

with 7i > and 72 > 0. We will next explain what restrictions are imposed on the functions Gi, G2, and 
F. The function G2 is analytic in the closed half-plane C+ = {A : i?e(A) > 0} and such that G2(0) ^ 0. It 
is also assumed that the function Gi is analytic in C+ except for a simple pole at A = with residue 1, that 
is, Gi(A) = J + G(A) where G is an analytic function in C+. In addition, we suppose that Gi is a nowhere 
vanishing function in C+ = {A : i?e(A) > 0}. Similarly, the function F is analytic in C+ and has a simple 
pole at A = with residue a > 0, that is, F{X) — j + F{X), where F is an analytic function in C+. We also 
suppose that 

Gi(A)^^G2(A)e^(^)| <exp{-|A|*} (16) 

as |A| 00 in C+, for some S > (much less is needed about the decay of the fmiction in (HH])). 

Remark 5.1 We assume that the function Gi does not vanish in C+ in order to justify the existence of 
the power Gj^ . Here we use the fact that for a nowhere vanishing analytic function / in C+, there exists an 
analytic function g on C+ such that /(A) = e^''''^ A G C+ (see Theorem 6.2 in [20]). 

The next result provides an asymptotic formula for the inverse Laplace transform of the function /. 

Theorem 5.2 Suppose that the functions Gi, G2, and F are such as above. Then the following asymptotic 
formula holds: 

M{y) = -l=«i+^G2(0)e^(")y-t+^e2v^v^ (l + O (j/-^)) (17) 

as y 00. 

Proof. Using the Laplace transform inversion formula, we see that for every e > we have 

^'^^y) = Iz=e+zr^i^)^^''^^- I* follows that M (ay) = L^e+tr ^i^)^"^''^^- By Cauchy's 
formula, we can deform the contour of integration into a new contour rj consisting of the following three 
parts: the half-line {—ooi,—y~^i]^ the half-circle F in the right half-plane of radius y~^ centered at (it is 



6 



oriented counterclockwise), and finally the half-line [y ^ijOoi). It follows that 



1 /■ ~ 1 /-oo 

h{y) + I2{y) + h{y)■ 



m 



We will first estimate l2{y)- This will give the main contribution to the asymptotics. By making a 
substitution z = ?/~2e**, < < ^, we see that 



I _ 1+21 f ^ 



exp i^F (y"5e»e^ | exp {a^i^e-*^} exp {aVye''^} e^^d6'. 
Next, taking into account the formula ^/y (e'^ + e^*^) = 2y/ycos9 ~ 2y/y + 2y/y(cos 6* — I), we obtain 

I,iy) = i-y-^e^-v^ /"^ e'^(i+-^)Gi {y-h:^'y G2 {y-h^') 
exp (y-5e»«^| exp{2a7^(cos6l - I)}d6'. 

It is easy to see that 

G2 {y-h'') - G2(0) = O (y-^) , exp {f (y-ie*«) } - e^W = O (y-^ 
on the contour F. Moreover, using (|43p and the mean value theorem, we obtain 

on r. Therefore, 

G, {y-h^'y'-y^^0{y^\e\)+0{y^) 
on r. It follows from (HH), and (HH) that 

Hy) - ^G2(0)e^(")y-^'^e2"v^ exp {2aVy(cos - 1)} 



(19) 

(20) 
(21) 



(22) 



1 + +oi\e\)'^d9. 



(23) 



We will next employ Laplace's method to estimate the integral appearing in Consider the integral 

g-'>'^i^)^(^x)dx, where $ G G°°[a, 6] and ?/; G G°°[a, 6] (much less is needed from the functions $ and ?/;), 

and assume that there is an a;o E {a,b) such that $' (xq) = 0, and <& (xq) > throughout [a,b]. Then the 
following assertion holds: 

Theorem 5.3 Under the above assumptions, with s > and s — > 00, 



A + o(l 



(24) 



where A = V^'tp (xq) {xo)) . 



The proof of Theorem 15. 31 can be found, e.g., in [20] . 

Using ([M)) with a = — ^, 6 = f, $(a;) = 1 — cosx, ip{x) = 1, xq = 0, and s = 2ay/y, we see that 



g2aVy(cose-l)^^^ V|y-1 ^ y^cx). (25) 



Similarly 

|5»|e2"v^(™-«-i)d6i = o(y-^) , y->oo. (26) 

Therefore, (gSl), (123, and (HH) give 

/2(y) = ^-^G2(0)e^(«)2/-t+^e2"v/s(l + o(y-i)), y ^ oo. (27) 
Moreover, using (fT6|) we obtain 

/i(2/)+/3(2/) = 0(exp{-c/}), y^oo, (28) 
for some c> 0. It follows from ((THl), dSZ]), and (gS]) that 

M(a2/) = ^-l=G2(0)e^(°)y'^+^e2"v^(l + o(y-3)) , y ^ oo. (29) 



Now it is not hard to see that ([29]) implies p?)) . 
This completes the proof of Theorem 15.21 



6 Proof of Theorem 12.41 

We will first discuss the properties of the following complex function: 

$s(z) = z cosz + s sinz , z G C, (30) 

where s > — 1. The next lemma concerns the zeros of $s. A special case of this lemma was stated in [12] 
without proof and was used in [7] and [19j . 

Lemma 6.1 For all s > —1, the function $s has only real zeros. 
Proof of Lemma \6.1i For every n > 1, put Pn{z) = z 

riLi i^-w)- The function P„ is a polynomial 



of degree 2n + 1, all of whose roots (z ^ k, k £ Z, \k\ < n) are real. Put Qn{z) = z "^^^^ {z'^Pn{z)) = 
sPn{z) + zPl^{z). Then Q„ is a polynomial of degree 2n + 1 which vanishes at z = 0. Also, by RoUe's 

theorem, the function 4- (•2*-Pn('S)) vanishes at points strictly between k and fc + 1, — n < k < n, since the 

dz 

function z''Pn{z) vanishes at those points. It follows from the previous considerations that Q„(z) has all its 

sm IT z s 
2n + l roots that are real. But P„(z) — > by the product formula. Hence, (3„(z) — > — sinvrz + zcosttz, 

TT ^ TT 

and the desired conclusion that all the roots are real follows from the Rouche-Hurwitz theorem. The proof 
above implicitly used the condition s > — 1 (for otherwise z*P„(z) does not vanish at the origin). The result 
for s = —I can be derived from that for s > — 1 by a limiting argument. 
This completes the proof of Lemma 16.11 

In the present paper, we consider only the case where s > 0. It is clear that the function $s is odd and 
satisfies $s(0) = 0. 

Definition 6.2 For s > 0, the smallest positive zero of the function will be denoted by r^. 
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The number plays an important role throughout the paper. It is not hard to see that tq = ^, and | tt 
as s — > oo. Moreover, the function s i— > is differentiable and increasing on (0,oo). Indeed, the value of 
for < s < oo is equal to the first coordinate of the point in where the segment y = —s~^x, ^ < x < tt, 
intersects the curve y = tanx. In addition, we have = 0^^(s), < s < oo, where (f){u) — — u(tanu)^^, 
^ < u < TT. It is also clear that sin(ro) = 1, cos(ro) — 0, and $q (tq) = ~f ■ Moreover, if s > 0, then 

sm{rs) > 0, cos(rs) < 0, and ^'^ (r^) < 0. (31) 

By Lemma l6.H the function ps defined by 

Ps{z) = z cosh z + s sinh z = —i^siiz) (32) 

has only imaginary zeros. 

The next lemma concerns the mixing distribution densities. 

Lemma 6.3 Suppose that a > 0, 6 < 0, and c > m the Heston model. Then the probability distribution of 
the random variable at in ^ admits a density mt for every t > 0. 

Proof. We will first prove that the probability distribution of the random variable Ygds admits a density 
mt. This will allow us to prove the existence of the mixing distribution mt, since mt can be determined from 
the formula 

rntiy) = ^mt (t-h^) . (33) 



It follows from Theorem 14.51 that 




A /■* . 11 f abt] I VWTX 

= cxp 



c 



V^2TAcosh(itVfo2 + A) - 6sinh(iiV6M~X) 



f j;oc-^Asinh(l^Vb^) 1 

^^^\ Vfo2TAcosh(i<V62 + A) - 6sinh(iiV62 + A) J ' 

Denote by and '^2 the functions on the right-hand side and the left-hand side of formula ([M)) . respectively, 
and put 

u,,, - -4i-V|,|,| (35) 

where t > and b < 0. It is not hard to see that the function can be continued analytically from the 
half-line to the half-plane 

Cb^t = {X e C : Re{X) > -b^ + ub,t} ■ (36) 

This can be shown using the properties of the zeros of the function ps defined by ([5^ and the fact that 
A = is not a singularity of the function ^'i. On the other hand, the function is analytic in the half-plane 
C+ = {A e C : Re{\) > 0} (use formula with A > to prove the finiteness of the derivative of the 
function in C_|.). It follows that formula holds for all A S C-h^t- 

Since i^ G Cb^t for all ^ € R^, equation ([M]) implies the following formula: 

e-*«"di^t(a;) = *i(zC), ^ £ R\ (37) 

where vt stands for the probability distribution of the random variable ^ Ysds. 
It is not hard to see that 

f z \^ i (C2z2-C3)sinhzl 
*i(A) - ci — — exp • 



Ps{z)J Ps{z) 



where z = ^t^/lP + \, s ~ ~b, ps is defined by (|32p . and ci, C2, and C3 are some positive constants. 
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Next, we observe that for z lying in a proper sector of the right-hand plane and z ^ cx), we have 



and 



This implies that 



, o , sinhz ^ ,^ , 

(C2Z2 - C3) = C2Z + 0(1) 

Psyz) 



= 2e-' + 0{\e-^'\) . 



Ps{z) 

|*i(A)| <Cicxp{-C2|A|5} 



(38) 



for A G Cb^t and A — > 00. 

It follows from ((38)l that the function ^ 1-^ |5'i(i^)| belongs to the space (K^). Taking into account 
P7p and using Plancherel's Theorem, we see that the measure vt is absolutely continuous with respect to the 
Lebesgue measure on R^. Therefore, the density fht exists. This implies the existence of the mixing density 
mt, and the proof of Lemma 16.31 is thus completed. 

We will next prove Theorem 12.41 The following lemma will be used in the proof. In this lemma, we use 
the notation introduced in (1301) and (l35l). 



Lemma 6.4 Let a > 0, 6 < 0, c > 0, and let Yt he a CIR process satisfying the equation 
dYf — {a + bYt) dt + c\/YtdZt and such that Yq — yo a.s. Then the following formula holds: 



~>-v 



y 2exp{(6^-U6,t)2/}mt f^^^j 



2t 



■ exp 



abt 



it^X + Ub,t 



exp ■ 



ityo (A + Ub^t - b"^) sinh {\t^\ + Ub^t) 



(39) 



The functions on the both sides of \39fl are analytic in the right half -plane Cq — {A G C : -Re (A) > 0}, and 
the equality in \39fl holds for all A G Co • 



Proof. It was established above that formula ([34|) holds for all A G Cb,t, where <Cb,t is the half-plane 
defined by (|36p . It is not hard to see that Lemma WM follows from (p4|) . Lemma [631 a-nd (I33p . 
We will next compute the residue Aq of the function 



A(A) = 



^A + wfc,* 



^A + Ub^t cosh {\t^\ + Ub^t) + \b\ sinh {\t^JX + Ub^t) 



ity/X + Ub^t 



at A = 0. Since 



{rim) 



(40) 



(41) 



it is not hard to see that 

Ao = lim AA(A) 



it|6| 



t|fc| 



1 + cos - rijifci sin (rit|b|) 



(42) 



In the proof of (gl]), we used and ([21 



10 



Our next goal is to apply Theorem 15. 21 to the Laplace transform in formula (|39p . The numbers 71 and 72 
and the functions Gi, G2, and F in Theorem l5.2l are chosen as follows: 71 = 0, 72 = and Gi(A) = ^A(A), 

where A and Aq are defined by (HOI) and g2D, G2(A) = a|^^ exp {-^}, and 



F(A) 



ityo (A + Ub^t - fe^) sinh (^i^A + Ub,t) 



Note that Gi is a nowhere vanishing function in C+.It follows from (|38p that the decay condition in (fT6|) is 
satisfied. The next lemma provides explicit formulas for the residue a of the function at A = and for 
the number F{0). 



Lemma 6.5 The following formulas hold: 



tyoVi 
2c2|pi| 



fe2+4t-VL| >0. 



F(0) = 



2c2p2 



(771 - (4<-V|,|,| + b^) m) PI + (4i-'r|,|,| + 6^) p2 



The constants rji, 772, pi, and p2 in and {44^ are given by 

cos (^it|fc|) 



P2 



128 



Pi 



it|b| 



771 = sm (^ri J , 772 
^2 



it|&| 



1 



cos (ri(|fc|) -ri 



t|b| 



sm r 1 



tt|b| 



1 + ^t\b\ j cos (r - r 1 ,|,| sin (r 1 ,|,|) + 2r , sin (r^ 



t|6| 



(43) 
(44) 

(45) 
(46) 

(47) 



Proof. We will need the first two coefficients in the power series representation for the function Ai^(A). 
It is not hard to see that 

(48) 



sinh ^-t^A + Ub^t^ = 1771 + 7?72A + 



where 771 and 772 are defined by (|45p . Next, using formula (j41[) and taking into account that 



it|6| l'it|6| 



— 2 sin 



{rit\b\) =-2sin(ri,|,|) 



we obtain 



*it|6| «oV^ + "&,t = PiA + P2A^ + 



(49) 



where pi and p2 are defined by and (|T7)l . Using (1^ . we see that ?7i > and 772 > 0. Moreover, 
(|49)) . and the definition of F imply the following equality: 



AF(A) 



t\b\ 



b^ + Xj (771+772A + ---) 



2c2 pi + P2A H 

tyo - 771 + {m {^t-'rltw + b') V2) A + • • • 



2c2 



Pi + P2A + 



(50) 



where pi, p2, rji, and 772 are defined in (|T7|) . and pS)) . 
Now, it is clear that Lemma 15751 follows from ([50]) . 
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We will next complete the proof of Theorem 12.41 and compute the constants appearing in it. By applying 
Theorem 15.21 to the Laplace transform in (1551) . we see that 

y"5 exp { {b^ - Ub^t) y] mt ^^^^y; a, b, c, yo^ 

«^-7^a|^^ exp |-^| e^(0)y-i+7^e2^V^ (l + O (y-^)) (51) 



2V^ 

as y ^ 00. Next, replacing y by y'^^ in formula (|5ip . we obtain 

4+ 



mt[y;a,b,c,yoj = Y'^a'^ ^ — — exp|--2-je 1^ 



_ 1 J_2a 

y ^T'T^'exp 



[V2^ 1 






I exp 1 







as y ^ 00. Now it is clear that ((52)) implies Theorem 12.41 In addition, the following lemma describes the 
constants appearing in Theorem 12.41 



Lemma 6.6 The constants A, B, and C in Theorem\2.4\ are given by 



B = c-W2at, and C = {2c^) t [b^ + 4t"^r1 

where the numbers Ao > 0, a > 0, and F{{)) are defined in \4-^ , and (^Tp, respectively. 

This completes the proof of Theorem 12.41 
Corollary 6.7 The following formula holds: 

mt {y, a, 0, c, yo) = Ae-^y" e^yy-i+^ (l + O (y-^^ 



asy^oo, where A = ''^^^jj ^ exp {^}, B 2v^^ and C ^ ^ 

Corollary 16 . 71 can be derived from Theorem l2.4l and Lemma 16.61 by using the equality tq = § and the fact 
that for 6 = we haye Aq = a = ^^"'I , and F(0) = — 

7 Proof of Theorem D 

We will next formulate a theorem concerning the asymptotic behavior of certain integral operators. This 
result is a minor modification of Theorem 4.5 established in [12j . 

Theorem 7.1 Let A, ^, andb be positive Borel functions on [0,oo), and suppose that the following conditions 
hold: 

1. The function A is integrable over any finite sub-interval of [0,oo). 

2. The function b is bounded and lim b(jj) — 0. 

y^oo 

3. There exist j/i > 0, c > 0, and 7 with < 7 < 1 such that Q and b are differentiable on [yi, 00), and in 
addition |C'(y)| < cy^^'Ciy) and \b'{y)\ < cy'^b{y) for all y>yi. 

4. For every a > 0, there exists ya > such that b{y)({y) > exp |— ay"'}, y > ya- 
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5. There exists a real number I such that A{y) — eJyC,{y){l + 0{b{y))) as y oo. 
Then, for every fixed fc > 0, 



J ^(2/)exp|- 



- + k^y^]}dy 



/TT 



The only difference between Theorem 17.11 formulated above and Theorem 4.5 in [12j is that in Theorem 
17.11 we do not assume the integrability of the function ( near zero. It is not hard to see that Theorem 17.11 
can be derived from Theorem 4.5 if we replace the function ({y) by the function yl(y)e~'^ near zero. 

Let A{y) — y^^mt{y)e^y^ , k = \Jc + I = B, ({y) ~ Ay^^'^'^ , and b{y) = y^^, where the constants 
A, B, and C are as in Theorem l2.4l Then, it is not hard to see that condition 2 in Theorem 1 7 . 1 1 follows from 
formula ([8]). In addition, it is clear that condition 3 with 7 = 1 holds. The next lemma shows that condition 
1 in Theorem 1 7 . 1 1 also holds. 



Lemma 7.2 For every s > 0, y ^mt{y)dy < oo. 

Jo 



Proof. The function on the right-hand side of (|39|) is integrable with respect to A over the interval 
(1, oo). Suppose h is any positive function on [0, oo) such that its Laplace transform has this property. Then 

Jo ^{y)y^^'^y < °° f^^' s > O. it follows from this fact and p9|) that the function y y^im,t ("^^v^) 
is integrable over any interval of the form [0, s] with s > 0. This implies Lemma [7jH 
Next, applying Theorem 17. II we see that 

^ y-^mt{y)e^^\^-{^ + ^-^y^dy 

= A^exp kl-^z-i+^ exp {s/fc-^V^} e-^^- (i + O (z-^)) (53) 

as z — > cx). Replacing z by in formula ([55)) and taking into account formula Q and the equality 

Vsc + t . 

k — ^ — , we obtam 

2\/2 

D.ixoe^^x) = -A_2^^^-^t-i^t.iSC + tri~^ 

{\ogx)-i+^ exp I ^1 ^g^^^^ 1 Vlog^l (l + o((logx)-3)) (54) 

as a; ^ 00. 

Now it is clear that formula (|54p implies Theorem 12.21 Moreover, the following lemma holds: 
Lemma 7.3 The constants Ai, A2, and A3 in Theorem \2.2\ are given by 

A 3, a 1 o , _ , „ r 52 

xoe- 



A, . _-,2-t+^r-if^(8C + r-^exp^^^^^ 



A2 = 1 — r, o-nd A3 ^ ^ + -i^l^il where A, B, and C are defined in Lemma\K 

ti(8C+t)i ^ ^vt 

Let p G M, and denote by Ip the moment of order p of the function Dt, that is, 
lp=E [Xf] = xPDt{x)dx. The next result was obtained in 1 . 

Lemma 7.4 For a > 0, b < 0, c > 0, and p € M., the following statement holds: < 00 if and only if 
2 — A3 <p<A3 — 1, where the constant A3 is such as in Theorem \2.2\ For b = Q, qp < 00 if and only if 



1 V47r2 + cH"^ 1 V47r2 + cH"^ 

< p < 



2 2ct 2 2ct 

It is not hard to see that Lemma [73] and more precise integrability theorems for the distribution of the stock 
price follow from Theorem [ 
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8 Proof of Theorem [23] 

Let Yf be the volatility process in model ([T|). Then is a squared Ornstcin-Uhlenbeck process. The Laplace 
transform of the law of the squared Ornstein-Uhlenbeck process was found by Wenocur and by Stein 
and Stein 'I95. Another explicit expression for this Laplace transform is given in the next formula: 



Vo 



exp 



exp ■ 



exp ■ 



exp 



A / Y^ds 



w 



j/pAsinh {ty/w) 



exp 

/w cosh {ty/w) + q sinh {t^/w) J 
■w?q^\ (sinh (ty/w) — ty/w cosh {ty/w)) 

W {y/w cosh {ty/w) + Q siuh (ty/w j) 

^g^A (4sinh2 (^^) - iViZ^sinh(tV^)^ 
wi {y/w cosh (ty/w) + gsinh (ty/w)) 



w cosh {ty/w) + q sinh {ty/w) ^ 

2niqyoX (cosh (tVw) ~ 1) 



y/w {y/w cosh (i^/w) + Q siii^ (^V^)) 



(55) 



where A > and w = q'^ + 2a^X. We will next sketch the proof of formula ([55]) . The first step is to replace 
the symbols S, 0, k, ctq, and A in formula (8) in [19] by q, m, a, yo, and tX, respectively, and take into account 
the following relations between the notation in [T^ and in the present paper: A = , B = ^ , C 

= , and ak^t 



= 2^y/fT2^, b = 



tV^' 



q_ p; _ rnq /-i_ A 

2 J- 2a^X. We also combine the terms ^f-^a^k^t and 



— i log 1^ ("f + 1) + ^ (1 ~ ) e^"'^ *| in the expression for N in formula (7) in [T3], and after somewhat 
long and tedious computations show that formula (8) in ^19j and formula ([55)) in the present paper are 
equivalent. 



Remark 8.1 Lemma [6.31 states that for the Heston model, there exists the mixing distribution density mt 
for every t > 0. The same statement holds for the Stein-Stein model. This can be established using formula 
()55p and reasoning as in the proof of Lemma [ 



It follows from Remark [8T] that 
'VO exp <j -A 



Y^dt 



'^^y mt{y)dy 



2Vt Jo 



'^yy ^^mt{t 2y2^ dy. 



(56) 



We will next obtain a sharp asymptotic formula for the mixing distribution density in the Stein-Stein model 
using formulas (|55|) . (|56l) . and the methods developed in Section [6] 

Recall that for s > 0, we denoted by the smallest strictly positive zero of the function ^s{z) — 

zcosz -t- ssinz. It is clear that zcoshz -I- ssinh^; — ~i^s{iz). For q > and t > 0, put Vq^t - 
([55)1 and dnD) give 



-0-. Now 



e-^yy- 



2 exp {(g^ - Vq^t) y} mt 



dy 



^/2^ 



■ exp 



where 



Fi(A) 



ity/X + Vq/t 
^qt {ity/X + Vq^t) 



Fi(A) exp {F2(A) + F3(A) + F4(A) + F5(A)} , 

1 

iylt (A + Vq^t - q^) sinh {ty/X + Vq^t) 



F2{X) 



2cr2$^j {it^X + Vq^t) 



Fz{X) 



imqyot (A -|- Vq^t - q^) [cosh {ty/X + Vq^t) - l] 



F4(A) 



a'^yJX + Vq,t^qt {ity/X + Vqj) 

ini^q^t (A + Vq^t - q^) [sinh {t^X + Vq,t) - t^X + Vq^t cosh [ty/X + Vq^t)] 

20-2 (A ^- Vq^t) ^qt {ity/X + Vq^t) 



(57) 

(58) 
(59) 
(60) 
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and 



^5 (A) 



Isinh^ — t-\/A + Vq^t sinh {t^\ + Vq,t) 



It is not hard to see that the functions Fi, F2, F3, F4, and F5 have removable singularities at 

^2 



(61) 



A = —V, 



q.t 



'_qt_ 



In addition, these functions are analytic in C+. Let us denote by Ai the residue of the function ^Fi at A = 0. 
It is not hard to see that 

Ai — „ 1 , r ■: r ; ^ rr . (^2) 

1(1 + qt) cos (r,t) - rqt sm (r,t)| 
Our next goal is to apply Theorem 15.21 to ((57l) . Put 



71 = 0, 72 = :^, Gi(A = ———===--, 

2 ^qt {lty/\ + Vq,t) 



G2(A) 



2i ^ i 

e 2 Aj , 



and 



F(A) = F2(A) + FsiX) + F4(A) + F5(A). 



(63) 
(64) 



In the sequel, the symbols aj and Fj (0) will stand for the numbers in the formulation of Theorem 
associated with the function Fj, 2 < j < 5. We will next compute these numbers. The following formula 
will be helpful in the computations: 



{vq^t -q^ + X)X _ {vq^t -q^)+X 



Tl + T2X + ■ 



where 



^qt (it^/X + Vq^t) Cl + C2A + • • ■ 

t"^ ((1 + qt) cos {rqt) - rqt sin (r,t)) 



(65) 



Cl 



C2 = 



2r, 



qt 



((1 + qt) cos [rqt) + rqt sin {rqt)) 



^rlt 



Tl 



2 [vq^t - q^) rqt 



t"^ ((1 + qt) cos {rqt) - rqt sin {rqt)) 



>0, 



and 



T2 



Cl - {vq,t - q^) C2 

Ar^t ((1 + qt) cos {rqt) - r,,* sin (r,t)) - [vq.t - g^) ((1 + qt) cos (r^t) + r,t sin (r^t)) 



Here we use the facts that v, 



q.t 



2rqtt'^ ((1 + qt) cos {rqt) - rqt sin {rqt)) 

-- < 0, COS {rqt) < 0, sin (r^J > 0, {rqt) = 0, 
(v) = (1 + 9*) cos (rqt) - r,t sin {rqt) , 



and 



*gt {Tqt) = -2 sin (v) - {rqt) = -2 sin (r,,*) . 



We will next employ ([551) - (pT|) to find explicit formulas for the numbers aj and Fj{Q) for 
2 < j < 5. It follows from and the formula 



sinh ( t^JX + Vq^t] = i sin {rqt) - i- — cos {rqt) X + 

zt, 



qt 
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that 



AF2(A) =_!^(ti+t2A + ---) ( i sin {rqt)-i^ cos {rqt)\ + 



= 2^Ti sin (rgt) + 7^ [t2 sin (r^t) - n cos (r^t) ) A 



2r, 



qt 



Therefore, 



and 



V t 

a2 = -^'^i sin {rqt) > 



^2(0) = |§Ti sin (r,0 + f r2 sin (r,*) - n cos (r,0 



Moreover, ((59l) . ([65]) . and the fact that 



2r. 



qt 



cosh (t^A + fg,f) - 1 _ .1 - cos (rgt) .t^ {2rqt sin (r^t) + cos (r^t) - 1) 



A + --. 



give 



AF3(A) 



20-2 



(n + T2A + ---) j 



1 - COS (rgt) .t^ {2rqt sin (r^t) + cos (r^t) - 1) 



m^q^t 1 — cos(r„() m^q^t 
-Ti ' 



2^2 



2a2 



1- cos (rgt) (2rgt sin (rgt) + cos (r„t) - 1) 

T2 Ti 



rgt 



2rg^ 



A + -- 

A + 



Hence, 



and 



as = 



q'^t 1- cos (rgt) 



2a2 



-Tl 



> 



rgt 



i^3(0) 



m'^q'^t r 1 - cos (rgt) (2rgt sin (rgt) + cos (rgt) - 1) 

T2 Tl 



2cr2 



rgt 



qt 



In addition, (l60|) . ([65]) . and the formulas rgt cos (rgt) + qt sin (rgt) = and 



sinh {t^/X + Vq^t) - ii/XT^t cosh (t^/X^^v^t) 

A + Wg,t 

t^{l + qt) sin (rgt) .f^ sin (rgt) (1 - gt - r2j 



A + --- 



imply that 



AF4(A) = ^^(ri+r2A + ...) 



t^{l + qt) sin (rgt) .t* sin (rgt) {l - qt - rl^) 



A + 



'gt 



It follows that 



and 



m^q^t^ [1 + qt)sm{rqt) ^ 



2a2 



'qt 



qH ( t'^sin{rqt){l-qt~rl^) ^[l + qt)sin[rqt) 
Tl 3 hr2 5 



2a2 



'qt 



'qt 
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Finally, we have 

4sinh2 _ t^X + Vq^t sinh {t^XT^ _ (r,t sin (r,*) - 4sin2 ^) 



' qt 



— ^ ( V sm (r^t) - 4sin — 1 + — — sm (r,t) - cos (r,t) 

■^'qt ^ ^ ' ■^'qt \'qt 



and hence 



where 



AF5(A)-^^(ri+r2A + ...) 



2(7 

qt 



'qt 



(rgt sin (rgt) 



4 sin 



2 '^gt 



it^ / 3 



2 ^ 2r^, 



■sin (rqt) - cos (r^t) 



Using the previous equalities, we obtain 



^5 



m'^q^t'^ 4 sin^ ^ - r,t sin (r^t) 

-Tl- 



2(72 



and 



m2q3i4 4 sin2 !ii _ gin (j-qt) 



2a2 



'9t 



-Tl 



3t^ 
2^ 



- (rqt sin (rqt) - 4 sin^ + TTs" f — sin (r^t) - cos {rqt)) 

qt ^ ' '^"qt \"qt J 



Our next goal is to prove that > 0. Indeed, we have 



4sin2 



— - rqt sin (rqt) = 2 sin ^ (2 sin ^ - r,t cos ^ 



(72) 

(73) 
(74) 

(75) 



Since j < ^ < ^ and x < tan x for |- < a; < ^ , inequality as > follows from ([7^ and ([75)1 . 

We can now apply Theorem 15.21 with the data given by and . Recall that we denoted by a the 
positive number given by 

a — a2 + as + a4 + (76) 

where Q!2, ^3, 04, and as are defined in (ISS)) . (ffO)) . and (|7^. respectively. We also denoted by F the 

function F2 + F3 + F4 + F5 (see dMl)), and put 



F(0) - ^2(0) + ^^3(0) + F4(0) + Fs(0) 



(77) 



where ^2(0), -^3(0), ^4(0), and Fs(0) are given by dHT]), dM]), and respectively. By applying 

Theorem 15.21 to (l57t . we see that 



y-ie.p{{q'-vq,t)y}mt (^^^^ = Xf e^^%'h'^^ {l + O (y-i)) (78) 

as ?/ ^ cx). Next, replacing y by y'^2^ formula ([75)) . we obtain 



TOt (y) = — e 2 Aj^e"^*^"^ exp 



27rcr 









> exp 1 


I ' "j 





(79) 



as ?/ ^ cx). Now, it is clear that formula ([75)1 implies Theorem l2.3l Moreover, the following lemma holds: 
Lemma 8.2 The constants E, F , and G in Theorem \2.3\ are given by 

^e*Afe^("), F = a-^V2^t, and G = {2a^yh {q^ + T 



E = 



27rcr 



'-'rlt) 



where the numbers Ai > 0, a > 0, and F{0) are defined in I16S\) . ^76^ 1, and jy?] ), respectively. 
This completes the proof of Theorem [ 
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9 Proof of Theorem [23] 

We will first show that it is possible to apply Theorem 17. II with A{y) = y^^mt{y)e'^y^ , k = y^G + |, I = F, 

Q{y) — Ey^^, and b{y) — y^^, where the constants E, F, and G are such as in Lemma [8.21 It is not hard 
to see that condition 2 in Theorem 17.11 follows from formula ([7]). In addition, it is clear that condition 3 
with 7 = 1 holds. The validity of condition 1 in Theorem 17.11 can be shown by reasoning as in the proof of 
Lemma 17.21 Here we use formula ([57|) instead of formula ([SQ]) . 
It follows from Theorem 17 . 1 1 that 



J y ^'^t(2/)exp|- ^^ + ^^|rfy 



(80) 



log X 

as z — > oo. Replacing z by — j= in formula (|80p and taking into account formula ([9]) and the equality 

V 2t 

VSG + t . 
fc = = — , we obtam 



2%/2 



i^\/2 



(logo:)-^ exp I ^1 Vlog^l ^ (l + O ((logx)-* j j (81) 

as X — > oo. 

Now, it is clear that formula ([5T|) implies Theorem 12. II In addition, the following lemma holds: 

Lemma 9.1 The constants Bi, B2, and B3 are given by Bi ~ xfei^-i- (8G + ^)~^ exp | 2(8G+t) 

1^2 = 1 ''^^ — r, O'nd -63 = 1+ , where the numbers E, F, and G are defined in Lemma \8.2\ 

tI(8G+t)5 2Vt - 

It is an interesting fact that Theorem 12.11 with m = is a special case of Theorem 12.21 This will 
be explained below. Recall that in model ((1]), the volatility process is the absolute value of an Ornstein- 
Uhlenbeck process. The following explicit representation is valid for the Ornstein-Uhlenbeck process Yt 
satisfying the stochastic differential equation dYt = q (^m — Yt^ dt + adZt: 



Yt (g, m, a, yo) = e-'^'yo + (l - e-«*) m + ae^'* / e'"dZ„ 

Jo 

(see, e.g., [TB], Proposition 3.8). Therefore, 

Yt (g, TO, a, yo) = Yt {q, 0, a, yo + (e«* - l) m) . (82) 

It is known that squared Ornstein-Uhlenbeck processes are related to CIR-processes. Indeed, it is not hard 
to see, using the Ito formula, that the squared Ornstein-Uhlenbeck process Tt — Yt {q,0,a,yo)^ satisfies 
the following stochastic differential equation: dTt = (cr^ — 2qTtdi) + 2a^/T^dZt. Therefore, the uniqueness 
implies that the process Yt {q,0,a, zq)'^ is indistinguishable from the CIR-process Yt (a'^,—2q,2a,Zo)- It 
follows from that 

Yt {q, TO, a, yof = Yt (a^, -2g, 2a, {yo + (e^* - l) to)') , (83) 

and hence in the case where to = 0, the mixing distribution densities corresponding to the processes on the 
both sides of ([55)1 coincide. Therefore, all the results concerning the distribution density of the stock price 
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process in model ^ can be reformulated for model ([T]) with rn = 0. For instance, formula (jlSp becomes 
exp{-Xy}y^^mt (t^'^y^;q,0,a,yo^ dy 

' qt\( ^Jq^ + 20-2A \ ' 



2%/te 



exp 



2 J \ + 2cr2Acosli(V'?^ + 20-2 A) + qwi\\{t^ + 2ct2A) 



exp • 



y2Asinli(V9^ + 2cr2A) 



^J(p + 2cr2Acosh(tv/g^ + 2(t2A) + 9sinh(fv/g^ + 2(t2A) 



Summarizing what was said above, we see that Theorem 12.11 with m = follows from Theorem 12.21 and 
formula 

The next statement concerns the moment explosion problem for the Stein-Stein model. 

Lemma 9.2 Let q > 0, rri > 0, cr > 0, and p G R. Then, the following statement is true for the moment Ip 
of the stock price distribution density Dt in model (QJ).- 

1 VscTt 1 VscTt 

L < oo <J=4> p < P < X H F ; 

^ 2 2^/^ 2 2Vi 

1 
2ct2 

It is not hard to see that Lemma [9^21 follows from Lemma [7l4l and formula ([IS 



^here C = — (ig^ + 



10 Proof of Theorem [33] 

Let K > 0, and let / and g be positive functions on the interval [K, oo). We will write f{x) ~ 5(2;), x — > 00, 
if there exist constants ci > 0, C2 > 0, and Kq > K such that the inequalities cig{x) < f{x) < C2g{x) hold 
for all X > Kq. The following notation will be used below: Vo{k) = Vq{K) (see ([TT|) for the definition of Vo). 
The next lemma was obtained in [13 . 

Lemma 10.1 Suppose that there exist positive increasing continuous functions ip and (p such that 
limfc-^oo = limfc_>oo 0(fc) = 00 and 

as K —t 00. Then the following asymptotic formula holds: 



as K ^ 00. 

For the model in ([2]), let -0 be a function such as in the formulation of Lemma [10. II Our next goal is to 
find a function (j) for which formula (f84|) holds. It follows from (fT2|) that 

Vo{k) = e--^ (xoe'-^)' [ / yDr (xoe^^y) dy ~ e'' H Dt (xoe'^^y) dy] . (85) 

Now it is not hard to see that ^ and ([55]) imply 

Vb(fc) » fc-3+^e^^^e-(^^-2)fc^ /c->oo. (86) 

Put 



m = Ji2A3-4:)k-2A2Vk+ Q- Jjlogfc-21ogV^(fc). 
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Then we have (j){k) w ^/k, and shows that condition (IM)) in Lemma [10.11 holds . Applying this lemma 
and the mean value theorem, we see that 



i{k) = ^ 

Vt 



l(A3-l)k-A2Vk+ Q-^^logfc~^(A3-2)fc- Q-^)logfc 

Next using the fact that \/l — h = 1 — + O (/i^) as /i | 0, we obtain p^ . The proof of ([H]) is similar. 
Here we use (O instead of ([6l). 
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